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The drag of massless fermions in graphene double-layer structures is investigated in a wide rage 
of temperatures and inter-layer separations. We show that the inhomogeneity of the dielectric back- 
ground in such graphene structures for experimentally relevant parameters results in a significant 
enhancement of the drag resistivity. At intermediate temperatures the dynamical screening via 
plasmon-mediated drag enhances the drag resistivity and results in an upturn in its behavior at 
large inter-layer separations. In a range of inter-layer separations, corresponding to the strong-to- 
weak crossover coupling of graphene layers, we find that the drag resistivity decreases approximately 
quadratically with the inter-layer spacing. This dependence weakens with a decrease of the inter- 
layer spacing while for larger separations we recover the cubic (quartic) dependence at intermediate 
(low) temperatures. 
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Introduction The discovery of graphene [T] [2], a 
monolayer lattice of carbon atoms, opened up possibili- 
ties for exploring new phenomena in fundamental physics 
[31 H] and for creating a new generation of electronic de- 
vice applications [SlIH]. Based on monolayer graphene, 
novel double-layer structures have been recently realized 
experimentally [THH] where massless fermions in two lay- 
ers are coupled only via many-body Coulomb interac- 
tion. Some current efforts in graphene physics focus on 
graphene double-layer structures (GDLS) [TUHTO] with 
the aim to find new electronic properties, which will 
emerge from the inter-layer many-body Coulomb interac- 
tion of massless, chiral fermions with a unique Dirac-like 
spectrum. 

Frictional drag [T7] between spatially separate electron 
layers provides one of the most powerful tools for the 
study of interaction effects. Recently, Coulomb drag in 
GDLS has attracted substantial theoretical fT51 - B5] and 
experimental [7^ attention. Despite much efforts, no gen- 
eral agreement has been reached between experimental 
observations and theory, and there is still no clear under- 
standing of the dependence of the drag on the inter-layer 
spacing and on the carrier density while the full treat- 
ment of the plasmon-mediated drag and of the dielectric 
inhomogeneity effect for realistic samples is still missing. 

In the present paper we calculate the drag resistiv- 
ity in GDLS in a wide rage of temperatures, T, and 
of inter-layer separations, d, using the finite-T polariz- 
ability and the finite- T nonlinear susceptibility for indi- 
vidual graphene layers. We focus on three main ques- 
tions. First, we investigate the effect of the dielectric 
inhomogeneity of the GDLS surrounding environment 
on the drag. Then, we calculate the contribution to 
the drag made by double-layer optical and acoustical 
plasmon modes and study the dependence of the drag 
rate on the inter-layer spacing. We show that Coulomb 
drag in GDLS immersed in a three layer nonhomoge- 



neous dielectric medium (see Fig. l(left)) is significantly 
larger than that calculated for the respective averaged 
homogenous background. This enhancement is observed 
for temperatures up to the Fermi temperature, Tp, and 
it becomes larger with increase of the inter-layer spac- 
ing. We find that at intermediate temperatures the dy- 
namical screening of the inter-layer Coulomb interaction 
results in the plasmon enhancement of drag, which is 
strongly pronounced at large inter-layer separations. Our 
calculations for d = 30 nm show an upturn in the drag re- 
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FIG. 1. (left) A graphene double-layer system immersed in 
a three layer dielectric medium. The solid lines with spacing 
d represent the active and passive graphene layers, 1 and 2, 
separating different materials with dielectric permittivities ei , 
£2, and £3. (right) The solid curve corresponds to the ratio 
of the double-layer screening functions e^{q) to ea{q), which 
are calculated within the static screening approximation, re- 
spectively, in GDLS with a homogeneous average dielectric 
permittivity £13 and in GDLS with a nonhomogeneous dielec- 
tric background, corresponding to the three layered medium 
of the left figure. The dotted, dashed and dot-dashed curves 
show, respectively, the intra- and inter-layer effective dielec- 
tric functions, eii(qd), £22(?rf), and £12 (gd) in units of £13. 
The electron density in each graphene layer is ni — n2 — 10^^ 
cm~^. 

sistivity at approximately 0.15Tp. We calculate the drag 
rate both in the regime of strongly {kpd < 1) and weakly 
{kpd^ 1) coupled graphene layers (fc^- is the Fermi wave 
vector). For inter-layer separations such that kpd ~ 1 
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the decrease of the drag resistivity with d is approxi- 
mately quadratic and it weakens with increase of d. In 
the weakly coupled regime we recover the strong d"* de- 
pendence of the drag resistivity, calculated at T = O.lTp 
within the static screening approximation. At T = 0.22^ 
the inclusion of the plasmon-mediated drag via the dy- 
namical screening weakens the drag dependence on the 
spacing and results in a d~^ behavior. 

Theoretical concept Frictional Coulomb drag in 
double-layer electron systems manifests itself when an 
electrical current with density Ji driven along the active 
layer induces, via momentum transfer due to inter-layer 
Coulomb interaction, an electric field E2 in the passive 
layer, which is an open circuit. The transresistivity, de- 
fined as pd = —E2/J1, is the direct measure of drag, 
which is studied in experiment. In terms of the diagonal 
intra-layer conductivities, cti and a2, and the off-diagonal 
inter-layer drag conductivity, cfd, the drag resistivity po 
can be obtained by inverting the conductivity tensor. As- 
suming, according to the experimental situation, that 
an ^ 171,2, we have pn ~ — o'd/o'iC2- The drag con- 
ductivity between two graphene monolayers can be cal- 
culated from the general drag formula, derived within 
the Boltzmann equation [26], the memory function [27] . 
and the Kubo [28l [29] formalisms in the lowest order of 
perturbation theory with small inter-layer interaction. In 
GDLS the drag conductivity is represented as 
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where A is the normalization area, and q are the trans- 
ferred energy and momentum from layer 1 to layer 2 
at temperature T . The dynamically screened Coulomb 
propagator 142(9,0) describes the charge density fluctu- 
ations that realize the electron-electron interaction be- 
tween the graphene layers. Within the random phase 
approximation ^12 (9, 0) can be obtained from a standard 
2x2 matrix Dyson equation as 



^12(9,0) 
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where the double-layer screening function is 

e(9,w) = ei{q,uj)e2{q,uj) - vi2{qfTli{q,uj)Wl{q,uj) (3) 

with the screening function and the Lindhard polar- 
ization function |30j of graphene monolayers given, re- 
spectively, by £1,2(9,'^) = 1 - «ii,22(9)n° 2(9,^^) and 
n]'_2(9,w). In the actual calculations of the drag resis- 
tivity we make use of the exact semi-analytical formu- 
las from Ref. [3T| for the finite temperature polarizabil- 
ity n]'_2(9, w) of graphene monolayers, including both the 
inter-chirality and the intra-chirality subband electronic 
transitions. The temperature dependent chemical po- 
tential p{T) is determined by the carrier density n from 



the equation Li2 (— exp(— /2(i)/t))— Li2 (— exp(/2(t)/i)) = 
n/2t^ where Li2 {x) is the dilogarithm function, p, = 
p/Tp and t = T/Tp. In the static screening approxi- 
mation we use the total Lindhard polarization function, 
Tl^{q, 0) — 2kp/TTv, in the static limit for small momenta 
q < 2kp. Here v is the velocity of the Dirac fermions. 

In double-layer structures the bare Coulomb inter- 
action is given by a 2 x 2 matrix, which, in general, 
represents three different electron-electron interactions, 
the intra-layer, 1111(9) and ^22(9), and the inter-layer, 
^12(9) = ^^21 (9), given by 



v^j{q,d) 



qh]{qd) 
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where i,j = 1,2 are the graphene layer indices. The 
effective dielectric functions, eij{qd), take into account 
the inhomogeneity of the dielectric background of the 
GDLS, which plays an important role in determining the 
bare Coulomb interactions. In the dielectric environment 
consisting of three contacting media with different, fre- 
quency independent dielectric permittivities £1, £2, and 
£3, the effective dielectric functions of GDLS are ob- 
tained from the exact solution of the Poisson equation 
for the Coulomb potential and are given by Eqs. (5)-(7) 
of Ref. [13] As seen in Fig. 1 (right) the dielectric inhomo- 
geneity modifies largely the behavior of the bare Coulomb 
interactions and, what is especially important, the be- 
havior of the screening function in momentum space. In 
the long wavelength limit all three interactions are deter- 
mined by the same effective dielectric constant, given by 
the arithmetic average of the top and bottom surround- 
ing media, £13 = (ei -I- £3)72, and does not depend on the 
dielectric constant, £2, of the inter-layer spacer. This is 
not true for the double-layer screening function e{q^ 0). It 
is seen in Fig. 1 (right) that the static screening of GDLS 
£s{q) = £(9,0) in the long wavelength limit differs essen- 
tially from the screening function, £^(9), calculated for 
GDLS, immersed in a homogeneous dielectric medium 
with an average permittivity £13. As we see below this 
effect results in an overall significant enhancement of the 
drag resistivity. 

Furthermore, the other important quantity in Eq. ([T]) 
is the quadratic response function of the charge density to 
an external potential, Ti[q,0). Due to the linear disper- 
sion of Dirac fermions, rj(9,0) is no longer proportional 
to the imaginary part of the individual layer polarizabil- 
ity as is the case in the weak scattering limit for the 
usual two-dimensional electron gas with parabolic disper- 
sion. Here in the actual calculations of the drag resistiv- 
ity we use the finite temperature nonlinear susceptibility 
for individual graphene layers from Ref. [53] where the 
carrier transport time Ttr is approximated by a constant. 
This approximation is well justified at low temperature 
T <t:Tp [131[21]. The energy dependence of the scatter- 
ing time can be important in the vicinity of the Dirac 
point when p{T) <C T [53]. Here we limit ourselves 
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FIG. 2. The efTect of the dielectric background inhomogene- 
ity on Coulomb drag in GDLS. The top pair of symbol sets 
shows the log-log plots of the drag resistivity as a function 
of scaled temperature in GDLS with d — 5 nm inter-layer 
spacing, immersed in a nonhomogeneous (the upper set) and 
homogeneous (the lower set) dielectric background with the 
parameters of GDLS of Fig. [l] The two other pairs of data 
sets in the mid and bottom of the figure correspond, respec- 
tively, to calculations for the d = 15 and 30 nm spacing. All 
curves have been calculated within the static screening ap- 
proximation. 
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FIG. 3. Plasmon enhancement of Coulomb drag of mass- 
less fermions in GDLS. The top and bottom bold curves rep- 
resent the log-log plots of the drag resistivity versus scaled 
temperature for the d — 5 and d — 30 nm inter-layer spac- 
ing, respectively. The solid curves correspond to calculations 
using the finite temperature exact polarization and nonlin- 
ear response functions while the dashed curves are calculated 
within the static screening approximation. The parameters 
are the same as in Fig. [l] for GDLS with a nonhomogeneous 
dielectric background. The thin curves correspond to the cal- 
culations for d = 5 nm in GDLS with a homogeneous dielectric 
background. 



to the important range of not very high temperatures, 
T < Tp. Therefore, the use of the constant scattering 
time approximation near the upper limit of temperatures 
T ~ Tp should be still justified. Here we assume also 
that the intra-layer conductivities are restricted by im- 
purity scattering and use cti 2 = e'^Ttr£i,2FVit)/^ where 
T]{t) = t dz\z\/ cos\L'{z+^{t)/2t) for the temperature 
dependent intra-layer conductivities. 

Results and discussions Based on the formalism de- 
scribed above we present here our numerical calculations 
of the drag resistivity, carried out in a wide range of 
inter-layer separations and temperatures up to the Fermi 
temperature Tp. In Fig. |2]we study the effect of the di- 
electric background inhomogeneity on drag of massless 
fermions by comparing the drag resistivity calculated for 
GDLS immersed in a nonhomogeneous dielectric back- 
ground with that calculated for GDLS in a homogeneous 
background. The presented results for three different val- 
ues of the inter-layer spacing, d = 5, 15, and 30 nm, show 
that the effect of the dielectric inhomogeneity is impor- 
tant in GDLS with realistic dielectric parameters, cor- 
responding to the experimental samples of Ref. [T] This 
new effect comes from the momentum dispersion of the 
effective dielectric functions for the three layer nonho- 
mogeneous dielectric medium (see Ref. [T3|. It changes 
the respective bare Coulomb interactions Q and reduces 
the screening function ( 3| in a nonhomogeneous dielectric 
background (c/. Fig. IT I that in its turn enhances the ef- 
fective inter-layer Coulomb interaction given by Eq. ([2]). 
We find that for d = 5 nm the drag resistivity in GDLS 
with a nonhomogenous dielectric background is larger by 



a factor of 3 than that in a corresponding homogeneous 
dielectric medium. This increase strengthens with the 
inter-layer spacing and for d = 30 nm the difference is a 
factor of 4.7. As seen in Fig. [2] the overall qualitative 
temperature dependence of the drag resistivity is not af- 
fected when we include the effect of the dielectric inho- 
mogeneity or when we vary the inter-layer spacing. At 
low temperatures the drag resistivity increases quadrati- 
cally with T due to the thermal broadening of the Fermi 
surface and at temperatures approximately T = 0.35Tp 
it shows a maximum. At higher temperatures the chem- 
ical potential /i(T) decreases and becomes smaller than 
T, the electron gases in the two graphene layers behave 
as a Boltzmann gas and the drag resistivity decreases 
approximately as the fourth power of temperature. 

One can see in Fig. [3] that the dynamical screening 
changes the temperature dependence of the drag resistiv- 
ity, especially, at intermediate and higher temperatures. 
For d = 5 nm the static screening approximation provides 
an adequate description of drag at low temperatures up 
to Q.2Tp. At intermediate temperatures the double-layer 
optical and acoustical plasmon modes become thermally 
excited and the drag resistivity increases a little by the 
plasmon-mediated drag in comparison with that obtained 
in the static approximation. At even higher tempera- 
tures, T ^ Tp, the dynamical screening has the opposite 
effect, it increases with T due to the temperature depen- 
dence of the dynamical polarizability and therefore the 
decrease of the transresistivity with T near Tp becomes 
stronger (c/. the solid and dashed curves). For the larger 
inter-layer spacing of d = 30 nm the static screening ap- 
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FIG. 4. The log-log plot of the drag resistivity as a function 
of the inter-layer spacing d for T — 0.2Tp (the solid and 
dashed curves) and for T = O.lTp (the dot-dashed and dotted 
curves). The solid curve corresponds to calculations using the 
finite temperature exact polarization and nonlinear response 
functions, the other curves are based on the static screening 
approximation. The dotted curve is calculated for GDLS with 
a homogeneous dielectric background, the other curves for 
GDLS with a nonhomogeneous dielectric background. The 
parameters are the same as in Fig. [l] 

proximation slightly overestimates the drag rate at low 
temperatures. One can see, however, that the dynami- 
cal screening causes an upturn in the drag resistivity at 
approximately = O.lSTp. With an increase of the 
inter-layer spacing the drag mediated by the electron- 
hole fluctuations decreases with d much faster than the 
plasmon-mediated drag (c/. the solid and dashed curves 
in Fig. |4| . Therefore the plasmon enhancement is more 
strongly pronounced for d ^ 30 nm than for rf = 5 nm. 
Note that the upturn temperature Tc in GDLS is smaller 
than that obtained for the usual two-dimensional electron 
gases [52] , 

In Fig. |4] we study the inter-layer spacing dependence 
of the drag resistivity at T = O.lTp and T = 0.2Tp. This 
dependence is mainly determined by the parameter kpd^ 
which is a measure of the inter-layer coupling. For small 
values of kpd 1 (for the density n considered here we 
have kpd ^ 1 for d ^ 6) the typical momenta q in drag 
scattering events are of the order of qt = T/v and do not 
depend on d. In this case we find that the drag resistivity 
decreases as pn d~^ with S <2. For larger separations 
with kpd ^ 1 the drag mediated by the electron- hole 
fluctuations is dominated by the momenta q ^ d~^ and 
we flnd that the drag resistivity calculated within the 
static screening approximation behaves approximately as 
[2S1I33]- At large separations the effect of the dielec- 
tric inhomogeneity on the spacing dependence of drag 
is weak ( cf. the dotted and dotdashed curves in Fig. |4] 
have almost the same dependence on d). However, the 
dynamical screening due to the double-layer plasmons 
weakens the inter-layer spacing dependence of the drag 
resistivity and we find that it behaves (the solid curve) 



approximately as d . 

In conclusion, we investigated the drag of massless 
fermions in GDLS in a wide range of temperatures and 
inter-layer separations. Our theory includes the effect of 
dielectric background inhomogeneity, which for realistic 
parameters results in a significant increase of the drag 
rate. At intermediate temperatures the thermally ex- 
cited double-layer plasmons cause an upturn in the drag 
resistivity. We find that the drag resistivity decreases 
with the inter-layer spacing approximately quadratically 
for inter-layer separations, corresponding to the strong- 
to-weak inter-layer coupling crossover. This dependence 
increases (decreases) with an increase (decrease) of the 
spacing. 
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